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Abstract 



We study n-point functions at finite temperature in the closed time path 
formalism. With the help of two basic column vectors and their dual partners 
we derive a compact decomposition of the time-ordered n-point functions with 
2" components in terms of 2" _1 — 1 independent retarded/advanced n-point 
functions. This representation greatly simplifies calculations in the real-time 
formalism. 
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I. INTRODUCTION 



Finite temperature field theory has been extensively studied for a large number of prob- 
lems in particle physics, solid state physics and the physics of the early universe. There 
exist real-time (RTF) and imaginary-time (ITF) formulations of finite temperature field 
theory The real-time formalism again comes in several variations [H. To study the 

relations between those different formulations of the theory has been an interesting sub- 
ject over the years [B— ITT[| . In Ref. [Q Evans gave relations between the RTF 3-point Green 
functions and their analytically continued ITF counterparts. In Ref. he noted that the sit- 
uation was more complicated for 4-point and higher order n-point functions. Taylor studied 
the corresponding relations for 4-point functions in ||12j| . In |i"3| , |l4|1 a column vector calculus 
for thermo-field dynamics was presented which simplifies the calculation of Feynman dia- 



grams for matrix valued propagators. More recently in Ref. |T5[ the relation between the 
time-ordered and retarded/advanced [|Il|] 3-point functions was reexamined in the closed 
time-path (CTP) formalism |3|,|T4||, and a simple decomposition of the 8-component real- 
time vertex tensor in terms of retarded/advanced 3-point functions was derived using outer 
products of the 2-component column vectors introduced in |13,|14|. It was shown in |T5|JT6[| 
that, due to orthogonality relations between the column vectors, this representation greatly 
simplifies calculations in the real-time formalism. 

The purpose of this short report is to establish relations between the time-ordered and 
retarded/advanced n-point functions by using the column vector technique. We will gen- 
eralize the work in Ref. |15| on the 3-point functions to n-point functions. In doing so we 
shall show that all 2 n time-ordered functions can be expressed through linear combinations 
of 2 n ~ 1 — 1 independent retarded/advanced functions. 

In Sec. |I] we will review the column vector representation of the two-point function and 



study how it is decomposed into retarded/advanced propagators. In Sec. |TTTj we will study 
the analogous decomposition for the 3-point functions. General n-point functions will be 
studied in Sec. |V|. Some conclusions are presented in Sec. [V]. 
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II. 2-POINT FUNCTIONS 



To establish our notation and for later use we first consider the single-particle propagator 
for a bosonic field theory. In real time, the propagator has 2 2 = 4 components since each of 
the two fields can take values on either branch of the real-time contour [|],|3|]. We generally 
follow the notation in Ref. [H|. 



The 4 components of the propagator can be combined into a 2 x 2 matrix 

(D n D l2 \ 

D=\ , (1) 

which can be rewritten 1 13|,14] in terms of an outer product of 2-component column vectors: 



Here 



n (Po) = a T\ ' n (~Po) = -(1 + n(p )) (3) 



is the thermal Bose-Einstein distribution, and Dr^a are the retarded and advanced propa- 
gators 

D R = D u - D 12 D A = D u - D 21 . (4) 
Their spectral representations are 

J-oo Itx Pq — io + le 

p~(p) is the (real) spectral density in terms of which all propagator components can be 
expressed via spectral integrals. 

Defining the following two column vectors 

en(p)=( 1 t^ ) ), e A (p) = -()), (6) 



n(p ) 



Eq. (H) can be written as 
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D{p) = -(D R (p) e A (-p)®e R {p) + D A (p) e R {-p)®e A {pj) . (7) 
We also define the "dual" vectors 

which satisfy (r 3 is the third Pauli matrix) 

e R ,A(p) = -r a e AR (-p) . (9) 

As a rule of thumb the vectors e RA are associated with legs which carry outflowing momenta 
while the e RA are associated with inflowing momenta. If the propagator D{p) = G^ 2 \p, —p) 
is represented by a Feynman diagram in which the momentum p flows from left to right, 
the left leg carries the inflowing momentum p (i.e. the outflowing momentum —p) while 
the right leg carries the inflowing momentum — p (i.e. the outflowing momentum p). This 
structure is reflected by the signs of the momentum arguments of the column vectors e R;A 
in Eq. ©. 



With the column vector contraction rule P~3|.|T"il 



(:)-C;)= oA+ °* < io > 

the basis vectors @ and their dual partners (||) satisfy the orthogonality relations 

e<*(p) ■ epip) = <W P = R,A). (11) 

Throughout this paper we will use latin letters a, b, c, . . . = 1, 2 for the usual thermal indices 
and greek letters a, /3, 7, . . . = R, A for retarded/advanced indices. 

The self-energy II, defined through the Schwinger-Dyson equation D~ x = Dq 1 — IT, is 
given by |L5| 



°b)-^)(>C^)-^«(-i^W-i)- (12) 



which can be rewritten as 
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n(p, -p) = -(n B (p) e fl (p)<g>e A (-p) + n A (p) e A (p)®e B (-p)) . (13) 

The two momentum arguments denote the two inflowing momenta in the legs connected to 
the self-energy. Let us define its R/A components via 

IW(p, -p) = (e a (p)®e p {-p)) • n(p, -p) {a, (3 = R,A). (14) 

The • denotes the contraction between the 2 vectors e in the outer product shown in Eq. (|H"D 
with the corresponding 2 dual vectors e carrying the same momentum argument in Eq. (|13|) . 
Using flTT| ) one finds 

Krr = Raa = , Il RA = U* AR = -U R = -U* A . (15) 



Note that this differs from the notation used in Ref. |TTJ who give Ura = ^-ar = and 
= Ur = H* AA . The reason for this discrepancy is that we use the convention that both 
momenta in H(p, —p) are inflowing, and our R and A indices indicate retarded and advanced 
boundary conditions for these momenta (see below). Aurenche and Becherrawy, on the other 
hand, use for the two-point function (contrary to their convention for the 3-point function 
which agrees with ours) one in- and one outgoing momentum. Following @ this implies 
that on the leg with the outgoing momentum the R and A indices should be reversed relative 
to our notation, in agreement with the above findings. 



III. 3-POINT VERTEX FUNCTIONS 

The RTF n-point Green functions have 2 n components G aia2 ^ an (pi,P2, ■ ■ ■ , Pn) where = 
1 or 2. We denote the truncated n-point function, from which the n external propagators 
have been removed, by r aia2 ... an (pi,P2, ■ ■ ■ ,Pn) and call it the n-point vertex. We adopt the 
convention that all momenta flow into the n-point vertex such that pi + P2 + . . . + p n — 
due to energy-momentum conservation. 

We will show here that the n-point vertex functions can be conveniently expressed as 
sums of outer products of n column vectors in a very similar way as Eq. (|13|) . For the 
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connected 3-point function such a representation was given in Eq. (35) of Ref. [15]. The 
analogous representation of the truncated 3-point vertex is obtained by reversing the signs 
of all lower components in the column vectors ||17|| : 



rOi,p 2 ,Ps) = (r abc (p 1 ,p2,P3) 



(1 + n 2 )(l + n 3 ) - n 2 n 3 



(1 + nO 
+ rR(pi,p 2 ,p 3 ) 



m 
m 

-1 -ni 



n 3 \ 
l-n 3 ) 

( ni 

1\ { n 3 \ 
1/ V-l-rJ 



F* Ri (pi,P2,P3 



(l + n 1 )(l + n 3 )-n 1 n 3 , . /l 

77— — v r B (pi,p2,p 3 ) 

(l + n 2 )-ri2 V-l, 



^2 
•1 - n 2j 



ni 

-1 -Tii 



"-2 
1 - n 2 



® 



1 



(l + m)(l + n 2 ) -nira 2 , v 



' 1 ' 



(l + n 3 )-n 3 ^ 
Here = n(p®), i = 1, 2, 3. Using (^) this can be rewritten as 

r(pi,P2,Ps) = rjt4A(pi,P2,P3) e R (pi)®e A (p 2 )®e A {p 3 ) 
+ Farr{pi,P2,P3) e A {pi)®e R (p 2 )(S)e R (p 3 ) 
+ Fara(pi,P2,P3) e A {px)®e. R {p 2 )®e A {p 3 ) 
+ Frar{pi,P2,P3) e R (pi)^>e A (p 2 )(S>e R (p 3 ) 
+ Faar(pi,P2,P3) e A {p 1 )®e A {p 2 )®e R {p 3 ) 
+ Frra{pi,P2,P3) e R (pi)®e R (p 2 )<g>e A (p 3 ) 

where 



n 3 

- n 3y 



(16) 



(17) 



r a /3 7 (pi,P2,P3) = (e a (pi)®e /3 (p 2 )®e 7 (p 3 )) •r(pi,p 2 ,p 3 ) (a, /3, 7 = -R, A) . (18) 

The indices R, A of the 3-point functions on the r.h.s. of (|PTD indicate whether a retarded 
or advanced propagator (D R or D A , respectively) is attached to the corresponding leg. 

As shown in Ref. ||11|| , the retarded/advanced 3-point functions on the r.h.s. of (|PT|) 
are straightforward analytical continuations iuj — >• p° + itj (J = 1,2,3) of the ITF 3-point 
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function, with tj > for each momentum corresponding to a leg labelled by R and tj < 
for each momentum corresponding to a leg labelled by A. Momentum conservation requires 
ei + ^2 + £3 = which implies Taaa = ^rrr = [[LJ]]. This explains why these two R/A 
functions do not appear on the r.h.s. of (|I7|). While Yaa...a — follows directly from the 
general identity || 

E r aia2 ... an = o, (19) 

a\,a,2,...,a n =l 



the vanishing of the n-point vertex with only retarded legs, jj = 0, is a result [Tj| of 
the KMS condition and is therefore only true in global thermal equilibrium. 
Comparison of ([17]) with ([16]) yields the identifications 

Tffi = Fraa, Ti? = Vara, F Ro = Taar (20) 

as well as the identities 

r , \ (l + ra 2 )(l + 7i 3 ) -ra 2 w 3 

rAHH(Pi,P2,ps) = 77- — x ^raa{puP2,P3) , (2ia) 

(l + m)-n! 
r /- \ (l + ni)(l + n 3 ) -nin 3 

rRAfl(Pl>P2,P3) = 77-, x T ARa(P1iP2,P3) , (21b) 

(l + n 2 )-n 2 
r / \ (1 +ni)(l + n 2 ) -nin 2 

rRRA(pi,P2,P3) = 77" x T AAr{P1iP2,P3) ■ (21c) 

(1 + n 3 ) - n 3 

The latter are manifestations of the general relation |18| 

n (1+nO 

r Q1 ... a „bi, . . . ,pn) = (-i)" +l 0K;= !! - a^dh, • • • ,Pn) 

{i|a,=A} 

n (i+rii)- n n< 
_l ij n (i+%)- n n/ Si -^ lPl '--- ,Pnj 

{i|»;=A} {i|c i= yl} 

which results from the KMS condition in thermal equilibrium. Here at = A, R if a = R, A, 
respectively. The second equality in ( p2|) is easily proven using the identity 
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IV. n-POINT VERTEX FUNCTIONS FOR n > 4 

Running through the manipulations of the previous Section in reverse order it is now 
easy to generalize Eq. (|16|) to higher order n-point functions. We illustrate the procedure 
for n=4; its generalization to arbitrary n is then straightforward. 

Using r AAAA (p 1 ,p 2 ,p 3 ,p i ) = = T RRRR (p 1 ,p 2 ,p 3 ,p i ), the analogue of (|T7|) reads 



= F RAA a(pi 


P2 


P3 


Pi) 


^r(Pi) 


»e A (p 2 ^ 


®e A (p 3 ) 


®e A {p±) 


+ ^arrr{pi 


P2 


P3 


Pa) 


£a(pi) 


^e R (p 2 ] 


®e R (p 3 ) 


®e R (p 4 ) 


+ T A raa(pi 


P2 


P3 


Pi) 


e A (pi)( 


»e R (p 2 ] 


®e A {p 3 ) 


®e A (p 4 ) 


+ F RAR r(pi 


P2 


P3 


Pa) 


£r(Pi) 


%e A (p 2 ) 


®e R (p 3 ) 


®e B (p 4 ) 


+ F A ara{Pi 


P2 


P3 


Pa) 


e A (pi)( 


®e A (p 2 ) 


®e R {j> 3 ) 




+ Frrar(pi 


P2 


P3 


Pa) 


ejj(pi) 


%>e R (p 2/ 


®e A {p 3 ) 


<8)e jR (p 4 ) 


+ F A aar{Pi 


P2 


P3 


Pa) 


e A (pi)( 


®e A (p 2 ) 


®e A {p 3 ) 


®Cfl(p 4 ) 


+ F RRRA (pi 


P2 


P3 


Pa) 


efl(pi) 


®e R (p 2 ) 


®e R (p 3 ) 


®e A (p 4 ) 


+ Trraa(pi 


P2 


P3 


Pa) 


6k(pi) 


2>e R {p 2 ) 


®e A {p 3 ) 


<8>e j4 (p 4 ) 


+ T A arr(pi 


P2 


P3 


Pa) 


(Pi) 


%e A (p 2 ) 


®e R (p 3 ) 


®eii(p4) 


+ F A rra{Pi 


P2 


P3 


Pa) 


(Pi) 


%e R (p 2 ) 


®e R (p 3 ) 


®e A (p 4 ) 


+ Fraar{Pi 


P2 


P3 


Pa) 


Cb(pi) 


%>e A (p 2 ) 


®e A (p 3 ) 


®e R (p 4 ) 


+ Tarar(pi 


P2 


P3 


Pa) 


6A (Pi) 


%>e R (p 2 ) 


®e A (p 3 ) 


®e R (p A ) 


+ ?rara(Pi 


P2 


P3 


Pa) 


efl(Pi) 


®e A (p 2 ) 


®e R (p 3 ) 


®e A {pA) 



The vertex functions satisfy the general definition 

r ai a 2 ...<*„(Pi,P2, •••,Pn) = (e Ql (pi)®e a2 (p 2 )®---®e an (p n )) • r(p 1 ,p 2 , . . . ,p n ) (24) 

(with «j = R, A) in terms of contractions with the basis vectors (||). Eq. ( p2|) gives the 
relations 



S 



r / \ (1 +n 2 )(l +n 3 )(l + n 4 ) - n 2 n 3 n 4 

rAflRH(pi,P2,P3,P4) = 7— s ^i? AAA (Pi iP2iP3iPa) , (25a) 

(l + m) - "i 

r / \ (1 + wi) (1 + w 3 ) (1 + w 4 ) - nin 3 n 4 

riMflfl(Pl>P2,P3,P4) = 77- s r Ai?AA(Pl,P2,P3,P4) , (25b) 

(1 + n 2 ) - n 2 

r / \ (1 +n 2 ) (1 + n 4 ) -nin 2 n 4 

riJ/Mfl(Pl,P2,P3,P4) = 77- v r AAKA(Pl,P2,P3,P4) , (25c) 

(1 + n 3 ) - n 3 

r / \ (l + ni)(H-n 2 )(l + 7i3)-ni7i 2 n 3 

r/?i?ifA(Pl,P2,P3,P4) = (1 +Wl ) _ Wl T AAAR (Pi, P2, P3, P4j , (25d) 

r / \ (l + 7i 3 )(l+n 4 ) -n 3 n 4 , , , , 

rAARfl(Pl,P2,P3,P4) = 7— 777- v r i j /UA (pi,p 2 ,P3,P4) , (25e) 

r / ^ ( X +^2)(l + n 3 ) -n 2 n 3 

^arra(Pi,P2,P3,P4) = 77- 777- l v ^ raarKPu P2, Ps, P4j , (25f) 

(1 + ni){l + n 4 ) - nin 4 

(1 + n 2 )(l + n 4 ) — n 2 n 4 

rAiUK(Pl ) P2,P3,P4) = 7— T7— v r RAi?A(Pl>P2,P3,P4) , (25g) 

(1 + m)(l + n 3 ) - nin 3 



which reduce the number of independent retarded/advanced 4-point functions to 2 
7. Combining (|23|) with (^) we get 

r(pi,P2,P3,P4) = 

TflAAA e fl (pi) ®e A (p 2 ) <&e A (p 3 ) <g>e A (p 4 ) 

(1 + n 2 )(l +n 3 )(l +n 4 ) - n 2 n 3 n 4 

+ t— v ^raaa e A (Pi ®e fl p 2 )®e fl (p 3 ®e fi p 4 

(l + ni)-ni 

+Ta_raa e A (pi ) <8) e/? (p 2 ) <8> e A (p 3 ) Oe^ (p 4 ) 

(l + m)(l+n 3 )(l+n 4 ) -nin 3 n 4 _ 

+ 7— v r AJMA e fl (pi)®e A (p 2 )(S)e fl (p3)(8)e fl (p4) 

l^i + n 2 j — n 2 

+r,4A/?A eA(pi)<8)eA(p 2 )®e fl (p 3 )(8)e^(p 4 ) 

(1 + n x )(l +n 2 )(l +n 4 ) -mn 2 n 4 

+ 7— r ^aara eR{Pi)®eR{P2)®e A (p 3 )®e R (p 4 ) 

(l + n 3 ) -7i3 

+r,4AAij e J 4(p 1 )<8)eA(p2)®e J 4(p 3 )(8)e R (p 4 ) 

, (l + ni)(l+n 2 )(l+n 3 ) -7iin 2 n 3 _ . , . , , 

+ 7— r T\ AAR e R (p 1 )®e R ($ 2 )®e R {p 3 )®e A {p A ) 

(1 + ni) -ni 

+r Mj4 A e R (pi)®e fi (p 2 )(g)e A (p 3 )(8)e A (j9 4 ) 

(1 + ra 3 )(l + n 4 ) — n 3 n 4 , n / \ ~ / \ / n 

+7TT — \7T~T — S r «i?AA Pi ®e A p 2 ^ (p 3 )®e R (p 4 ) 

(1 + ni){l + n 2 ) - n x n 2 

+~Tr AA r e R (pi)®e A (p 2 )®e A (p 3 ) ®e R {p±) 

(1 + n 2 )(l + n 3 ) — n 2 n 3T _ 
+ (1 + +n 4 ) - nin4 T nAAReA(p^e R (p 2 ^e R ( P3 )®e A ( P4 ) 

+^rara(pi,P2,P3,Pa) e R (p 1 )<g>e A (p 2 )®e R (p 3 )<g>e A (p4) 

9 



4-1 



+ nt n int n4 l ^ RAra &A(pi)9eR(p2)9B A (ps)9e B (p A ) . (26) 
(1 + n 1 )(l + n 3 ) -mn 3 

The generalization is obvious. We have chosen the R/A vertices with the smallest number 
of R indices as independent n-point functions; among them the functions with only a single 
R index are the fully retarded linear response functions ||18|| . Eq. (|22f) permits, however, to 
select any other set of 2 n_1 — 1 independent n-point functions and to rewrite Eq. (|2~6"D and 
its generalizations accordingly. 

V. CONCLUSIONS 



We have given relations for arbitrary n between the real-time thermal n-point functions 
in the single-time representation and the retarded/ advanced n-point vertex functions intro- 
duced by Aurenche and Becherrawy [|TTJ. These relations were expressed in a very compact 
way in terms of an orthogonal basis of two 2-component column vectors and their dual 
partners. The advantage of the R/A representation is that it diagonalizes the single particle 
propagator, shifting the thermal distribution functions from the propagators to the vertices 
In our column vector representation the thermal distribution functions appear only in- 



side the column vectors. Their orthogonality properties implement certain relations among 
the thermal distribution functions which lead to immediate cancellations between various 
different contributions |H| to complex Feynman diagrams; in other approaches to real-time 
finite temperature field theory these calculations often occur at a much later stage of the 
calculation. Thus our representation leads to appreciable simplifications for the Feynman 
diagram calculus in the real-time formulation of thermal field theory. 

Our column vector representation gives the thermal components of the time-ordered n- 
point functions as linear combinations of 2 n ~ 1 — 1 independent retarded/advanced n-point 
functions and their complex conjugates. Although we have here considered only bosonic 
n-point functions, the generalization of our results to fermionic fields is straightforward and 
involves only the simple substitution n —> —rif in the column vectors. The useful relation 
between thermal distribution functions given in the Appendix has been presented for the 
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general case, too. Similar relations as those given here for the truncated n-point vertex 
functions hold for the connected n-point Green functions; the only change is that in that 
case the signs of all lower components in the column vectors e must be reversed. 
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APPENDIX A: A USEFUL RELATION 



Let us write the thermal equilibrium distributions as 



n(xi) = — , n(-Xi) = -rji - nfa) . (Al) 

e x% — rji 



where Xi = (3p®, rji — 1 for bosons, and rji = — 1 for fermions. Defining 

^ 1 

t Xl) ~ ex P (xti Xi ) - nti n 



one can prove the identity 

na + vMxi)) - fi vMxi) = u ;=\ n{Xt \ . (a 3 ) 

i=i i=i ^12^1=1^^ 

It is easy to check the identity for k = 2 and prove it for k > 2 by induction. 

Due to fermion number conservation, an n-point function with n external legs has always 
an even number of fermionic legs. This implies 

fb* = l- ( A4 ) 
i=i 
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Since momentum conservation also requires 

n 

E^ = °> ( A5 ) 

8=1 

Eq. (TA^) yields the identity 

n n 

H{l + r H n(x i ))-'[[n(x i ) = (A6) 

i=i i=i 

if the product goes over all legs of an n-point function. Previously used identities like 
Eq. (34) in Ref. |j| are simple consequences of Eq. (|A6|). The physical interpretation of 
Eq. ( |A6| ) is that the probability for creating and for destroying n particles in the heat bath 
are equal. 
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